Direct, Loss- Tolerant Characterization of Nonclassical Photon Statistics 
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We experimentally investigate a method of directly characterizing the photon number distribution 
of nonclassical light beams that is tolerant to losses and makes use only of standard binary detectors. 
This is achieved in a single measurement by calibrating the detector using some small amount of 
prior information about the source. We demonstrate the technique on a freely propagating heralded 
two-photon number state created by conditional detection of a two-mode squeezed state generated 
by a parametric downconverter. 
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The photon number distribution is a key characteristic 
of every optical field. Several indicators of nonclassicality 
are based directly on the measurement of these statistics, 
including negativity of Mandel's Q Parameter and the 
negativity of the Glauber-Sudarshan P function . Fur- 
thermore many quantum information applications rely 
on light sources with well-defined photon number dis- 
tributions. For example, quantum cryptography based 
on single-photon protocols |2| requires the complete sup- 
pression of multi-photon components in order to guaran- 
tee security over longer distances. The ability to directly 
measure the photon number distributions, then, is impor- 
tant both fundamentally and technologically. Two main 
obstacles have hindered progress in such measurements: 
first, photon- number-resolving detectors were only de- 
vised recently |^|^; second, all available single-photon 
sensitive detectors exhibit finite efficiencies, such that in- 
trinsic losses often mask the signatures of nonclassical 
states. 

Currently there are three approaches to retrieving 
photon number distributions: using photon-number- 
resolving detectors via state reconstruction from 
homodyne tomography 0|; and using binary 

("click-counting'') detectors like avalanche photodiodes 
(APDs) [3, Q, la|- However, all of these approaches are 
compromised by loss and detector inefficiencies, which 
cause instabilities in the inversion algorithm used to re- 
construct photon number distributions from count statis- 
tics. This makes it difficult to reconstruct the quantum 
state of the source since the detector efficiency must be 
known accurately. Several photon-number-resolving de- 
tectors have high detector efficiencies, though these are 
usually accompanied by noise from dark counts, which 
also affects the fidelity of the inversion. In homodyne 
detection the calibration is made yet more difficult by 
the need to match the modes of the quantum state with 
that of the local oscillator. Two distinct approaches 
are currently known using APDs: measurements of the 
mean count rate as a function of beam attenuation 0, |^ 
and mode-multiplexing to implement a photon-number- 
resolving detector jsj . Photon number characteristics for 



classical states have been measured with APDs in the 
past 0j IS 0| , but the problem is much more intricate for 
nonclassical states, which do not maintain the form of 
their statistics when attenuated. Absolute APD quan- 
tum efficiency measurements usually require an indepen- 
dent measurement in which the detector response can be 
distinguished from external losses. Poor calibration may 
compromise the accuracy of methods such as the atten- 
uation approach that rely on well-known losses. In con- 
trast, the mode-multiplexing approach greatly increases 
the accuracy of detector calibration and does not require 
an independent measurement. 

In this Letter we propose and demonstrate a new APD- 
based approach to direct loss-tolerant photonic state 
characterization. We exploit partial a priori information 
about the photon number distribution, typically known 
from the state generation process, to accurately calibrate 
the total loss in the channel taken by the state of inter- 
est (hereafter referred to as the signal). Thus we gen- 
eralize the idea of a self-referencing detector, originally 
introduced in 1977 by Klyshko to different types of 
multi-photon states. We use a time-multiplexed detector 
(TMD) to emulate a photon-number-resolving detector. 
The measured count statistics p{k) are related to the pho- 
ton number distribution p{n) by p{k) = C • L(77) • p{n), 
where L(?7) is the matrix describing the binomial process 
of loss with an overall efficiency of rj and C is a matrix 
that takes into account that the TMD can only resolve 
up to a finite number of photons at the input [ll|||l8l|. 
Thus the photon number distribution at the source can 
be reconstructed from the count statistics by directly 
inverting these matrices or by using a maximum likeli- 
hood technique. We emphasize that we utilize this cali- 
bration technique to accomplish, from a single measure- 
ment set, both loss estimation and a reconstruction of 
the photon number distribution at the source - allowing 
loss-independent state characterization for a given gen- 
eration process without the need of a pre-measurement 
of the loss. As with Klyshko's original detector calibra- 
tion scheme our approach relies on knowledge that we 
infer from the state generation; previous theoretical and 
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FIG. 1: (Color online) The experimental setup: A mode- 
locked titanium sapphire laser (Ti:Sa) pumps sum frequency 
generation (SFG). The Ti:Sa is eliminated with Schott glass 
filters (not shown) and the SFG bandwidth is restricted with 
an interference filter (IF). This is used to pump PDC in a 
waveguide and the blue is removed with a dichroic mirror 
(DM) and Schott glass fihers (not shown). The PDC pho- 
tons are split at a polarizing beamsplitter (PBS). The trigger 
arm (refiected) contains a spectral filter (SF) and can be de- 
tected with either a single APD or a beamsplitter (BS) and 
two APDs detected in coincidence. The signal arm (transmit- 
ted) is analyzed with the TMD. The additional loss (ND) was 
optionally placed prior to the TMD. 



experimental work indicates that such inferences are rea- 
sonable_P,0, 13]. Wc emphasize that earlier theoretical 
work ,14| has shown that for known losses such compen- 
sation procedures are always possible when the overall 
detection efficiencies exceed 50% (and yet smaller effi- 
ciencies can be tolerated for specific classes of states). 

In general, any type of prior information can be used, 
but in this Letter we use the strict photon number corre- 
lations between the modes of a two-mode squeezed state; 
the state of the field is 



|V) = v/r~^^A"|n),|n), 



(1) 



where the modes are labeled by orthogonal polarizations, 
h and v, and A is the parametric gain Q]. The state we 
wish to characterize is the horizontally polarized mode 
(signal arm) that is conditionally prepared by detection 
of the vertical mode (trigger arm). This would normally 
require an independent measurement of the detector pa- 
rameters and channel loss, but using our technique can 
be inferred directly from the state characterization data. 
Historically, detector calibration was accomplished by as- 
suming only the first two terms of the sum in Eqn.^con- 
tribute, measuring each mode with an APD and compar- 
ing the singles counts to the coincidence counts. We ex- 
tend this idea in two ways: first, we incorporate k terms 
of the sum in Eqn. ^ where k denotes the number of 
photons detected as a trigger, allowing the use of higher 
parametric gains and the characterization of a broader 
range of states; second, we detect the signal with the 
TMD, allowing us to see the complete count statistics of 
the signal, from which we then derive the losses. 

The experimental setup is shown in Fig. ^ A mode- 



locked titanium sapphire (Ti:Sa) laser (100 fs pulses at 
800 nm and a repetition rate of 87 MHz) pumps second 
harmonic generation (SHG), which we subsequently filter 
down to a 2 nm bandwidth (FWHM). Type II PDC is 
generated in a 12 mm long z-cut KTP waveguide with 5 
/iW of second harmonic power. The orthogonally polar- 
ized daughter photons split into different spatial modes 
at a polarizing beamsplitter (PBS). The trigger arm is 
spectrally filtered with a 15 nm filter (FWHM). Details of 
this high brightness, waveguided downconversion source 
are presented elsewhere |l5[ . For single photon genera- 
tion, the trigger is a multimode fiber-coupled avalanche 
photodiode (APD) (Perkin-Elmer SPCM-AQR-13). Al- 
ternatively, for two photon state preparation, an addi- 
tional BS is inserted into the trigger arm and two APDs 
are detected in coincidence (3 ns coincidence window) as 
a trigger. A trigger is only accepted if it occurs in a well- 
defined 700 ps window relative to the Ti:Sa pump pulse. 
The signal arm is detected with the TMD. 

Our results are presented in Fig. [3 We begin with the 
case of single APD trigger (Fig. El a&c). Note that if 
an APD was used in the signal arm this setup would be 
equivalent to that of Klyshko. However, the use of the 
TMD allows us to verify the complete count statistics 
of the conditionally prepared state that can be used to 
verify the validity of of our assumption that A is small 
such that higher order terms of Eqn.Qare negligible. We 
define the "Klyshko efficiency" for our detector as[l9| 



Vk = p(chck|click) — 



where Ni is the number events where i photons would be 
registered by a TMD and p(click|click) is the probability 
of registering a click in the signal arm conditional on 
receiving a click in the trigger arm. Using the TMD, 
we are able to define the efficiency associated with single 
photon triggers as 



r?=p(l|t= 1) 



where p{i\t ~ 1) is the probability that i photons were 
registered in the signal arm given a single APD click in 
the trigger arm. This relation and the further relation, 
p{Q\t — 1) = 1 — can be used to deduce an overall signal 
efficiency of 37.3 ± 0.1% from the data. Accounting for 
losses, this corresponds to a single photon conditional 
preparation efficiency of 97%. 

In our approach we utilize the complete conditional 
statistics to ascertain 77 from the experimentally inde- 
pendent measurements of the different photon numbers 
n. Note that in the general case of k photons detected 
in the trigger arm we expect p{n < k\t — k) = Q if 
there is no loss or detector inefficiency (77 = 1), due to 
the prior information of number correlations in the two 
modes. Thus all probabilities p{n < k\t = k) > Q are 
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FIG. 2: (Color online) Count statistics detected (logarith- 
mic scale) using (a) single and (b) double APD as a trigger. 
The solid line (used to guide the eye) shows the Poisson dis- 
tribution with the same mean photon number as the data 
(0.376 and 0.623, respectively). This line illuminates the sub- 
Poissonian nature of our measured statistics. The photon 
number distribution obtained by using a maximum likelihood 
inversion with constraints p{n) > for (c) single and (d) dou- 
ble APD trigger by taking into account the efficiencies, which 
were 37.3% and 31.5%, respectively. 



caused solely by losses in the signal arm with 



p{n<k\t^k)= (^)?7"(l-?7) 



fc-r 



(2) 



In this way we can exploit all such contributions with 
n < k to obtain a value for the efficiency independent 
of all other experimental parameters like the parametric 
gain A and the loss in the trigger arm. We note that any 
state of light that has p{n) = for at least one value of 
n (such as a single mode squeezed state where p{n) — 
for all odd values of n) is a perfect candidate for this 
technique. 

For the double trigger (Fig. I^b&d), one can calculate 
the losses in the signal arm in three different ways us- 
ing Eqn. 121 Given a two photon detection in the trigger 
arm {t = 2) and a signal efficiency rj, it is straightfor- 
ward to calculate the following relations between these 
probabilities and the efhcicncy in the signal arm: 



(0) 



1 - ^p{0\t = 2) 



= -[l-v/l-2p(l|i = 2)] 



V 



(2) _ 



Vp(2|t = 2) 



(3) 
(4) 
(5) 



To evaluate the efficiencies tj'-'' from our raw data we 
must consider not only the losses in our case, but also 



the limitations of the photon-number-resolving capabil- 
ities of the TMD. This is done by multiplying our click 
statistics by the inverse of C and then applying Eqns. 3-5. 
We emphasize that this matrix can be obtained from the 
same measurement as the characterization and no sup- 
plemental measurement is necessary. Using the above re- 
lations we find that r]^°^ = 31.5±0.1%, rj^^^ = 31.0±0.2%, 
and 77^^^ = 32.1 ± 0.2%. This corresponds to an average 
efficiency of 31.5 ± 0.2% and a two-photon conditional 
state preparation efficiency of 90%. The reason that 
these numbers differ by a degree larger than the error 
is because the state is not a true two-photon state, but 
contains small contributions from higher photon num- 
bers. Qualitatively, it is simple to see that with higher 
photon-number contributions r;'-^-' would be larger than 
its true value because p{2\t = 2) would be higher than if 
p{n > 2) = 0. Simulations of our technique confirm this 
quantitatively. We note that this problem occurs in all 
calibration approaches using twin photon beams. How- 
ever, in our case we are able to assess directly the valid- 
ity of our assumptions by having access to the full count 
statistics, something that is impossible without photon- 
number resolution. Even with this consideration, our es- 
timate of the signal arm efficiency (mean of above num- 
bers with error bars encompassing the spread of values) 
is the most precise direct calibration available. 

To test the accuracy of our loss estimation a neutral 
density filter was inserted into the signal arm before the 
TMD, using the double- APD-trigger configuration. The 
filter was calibrated to have a transmission of 13.5±0.1% 
using the Ti:Sa laser and a linear photodiode. Using 
the previously mentioned loss relations a measurement 
of the efficiency was performed with and without the ad- 
ditional filter; The ratio of these two efficiencies gives a 
filter transmission of 13.8 ± 0.1%. The slight discrep- 
ancy in transmittance is likely due to the spectral dif- 
ferences between the laser and the PDC signal photons. 
(Non-degenerate PDC was used resulting in a spectrum 
with a different central wavelength and bandwidth from 
the laser. This calibration discrepancy reemphasizcs the 
need for a more accurate way to calibrate loss than with 
a prior measurement using classical light.) As is expected 
for this case of very high loss (95.5%) the fidelity of the 
inverted distribution decreases to 75%. Experimental at- 
tenuations performed with the single APD trigger were 
able to accurately reconstruct photon number distribu- 
tions (showing > 90% fidelity to a single photon state) 
for total signal loss of up to 88.0%. 

An important case to study is when the prior informa- 
tion used (e.g. photon number correlations) is not valid. 
The most extreme situation would be the substitution of 
uncorrelated light sources for the PDC source. We in- 
vestigate this effect theoretically by mixing in a coherent 
state into the PDC signal arm with various mean photon 
numbers. We find that the fidelity of our photon num- 
ber reconstruction stays above 90% for means up to 0.3 
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photons, independent of the loss; the additional mean 
number can be significantly higher than 1 photon if the 
efficiencies are higher than 55% . Another key obser- 
vation is that decreasing prior knowledge accuracy {i.e. 
increasing the mean photon number of the additional co- 
herent state) leads to inconsistent results among the ef- 
ficiency measurements. We confirm this experimentally 
for the extreme case of completely uncorrelated light by 
pumping a waveguide that is not periodically-poled and 
therefore produces no PDC but creating spurious fluo- 
rescence counts. This yields count statistics with a very 
low mean photon number, which we evaluate according 
to our previous analysis. Using the false assumption 
of photon number correlations results in: 1) inconsis- 
tent loss measurements and 2) unphysical reconstructed 
photon number distributions (distributions with negative 
probabilities) . 

Finally, to prove nonclassicality, we show that both 
our detected statistics and the inferred photon number 
probabilities result in negative values of Mandel's Q pa- 
rameter, which is a sign of nonclassicality. All Fock states 
result in Q = — 1 and all coherent states have Q = The 
count statistics detected using a single APD trigger give 
Q — —0.36 and the inferred probability distribution re- 
sult in Q — —0.97. Moreover, using two APD triggering, 
the detected statistics yield Q = —0.32 and the inverted 
distribution results in Q = —0.93. To test the consistency 
of this nonclassicality we also investigated the negativity 
of the P function for our data. This measure can be for- 
mulated in terms of conditions on the photon number 
distributions as 0,0] 

B{n) = {n + 2)p{n)p{n + 2) - {n + l)[p{n + 1)]^ < 0, 

where the inequality need only be satisfied for one value 
of n to assure the negativity of the P function. Our de- 
tected statistics do not satisfy the conditions, but our 
inferred photon probability distributions do. For the sin- 
gle [double] APD case, Si(0) = -0.13 [^2(0) = -0.11], 
once again showing the nonclassicality of the states. 

In summary, we have demonstrated a method of non- 
classical state characterization using mode-multiplexing 
and standard APDs. An inversion of the counting statis- 
tics was used in conjunction with a 'calibration-free' mea- 
surement of the loss in order to recreate the photon num- 
ber distribution of a heralded photon source based on the 
two mode squeezed state. In the future, this technique 
could be extended to more general nonclassical states 
where distinct properties which change under the influ- 
ence of loss are known. In the context of conditional state 
preparation our characterization technique can also be 
extended where loss calibration is not obvious. By utiliz- 
ing properties of known unconditioned statistics we can 
estimate the losses in the system and thus obtain the loss- 
tolerant calibration for post-selected conditioned subsets. 
Hence we expect that our detection scheme will become 



particularly relevant for quantum information protocols 
such as entanglement distillation 17]. 
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